Exact Solution of the Munoz-Eaton Model for Protein Folding

A transfer-matrix formalism is introduced to evaluate exactly the partition function of the Munoz-Eaton model, relating the folding kinetics of proteins of known structure to their thermodynamics and topology. This technique can be used for a generic protein, for any choice of the energy and entropy parameters, and in principle allows the model to be used as a first tool to characterize the dynamics of a protein of known native state and equilibrium population. Applications to a $\beta$-hairpin and to protein CI-2, with comparisons to previous results, are also shown.Comment: 4 pages, 5 figures, RevTeX 4. To be published in Phys. Rev. Let

Origin of entropy convergence in hydrophobic hydration and protein folding

An information theory model is used to construct a molecular explanation why hydrophobic solvation entropies measured in calorimetry of protein unfolding converge at a common temperature. The entropy convergence follows from the weak temperature dependence of occupancy fluctuations for molecular-scale volumes in water. The macroscopic expression of the contrasting entropic behavior between water and common organic solvents is the relative temperature insensitivity of the water isothermal compressibility. The information theory model provides a quantitative description of small molecule hydration and predicts a negative entropy at convergence. Interpretations of entropic contributions to protein folding should account for this result.Comment: Phys. Rev. Letts. (in press 1996), 3 pages, 3 figure

The correlation of RNase A enzymatic activity with the changes in the distance between Nepsilon2-His12 and N delta1-His119 upon addition of stabilizing and destabilizing salts.

The effect of stabilizing and destabilizing salts on the catalytic behavior of ribonuclease A (RNase A) was investigated at pH 7.5 and 25 degrees C, using spectrophotometric, viscometric and molecular dynamic methods. The changes in the distance between N(epsilon2) of His(12) and N(delta1) of His(119) at the catalytic center of RNase A upon the addition of sodium sulfate, sodium hydrogen sulfate and sodium thiocyanate were evaluated by molecular dynamic methods. The compactness and expansion in terms of Stokes radius of RNase A upon the addition of sulfate ions as kosmotropic salts, and thiocyanate ion as a chaotropic salt, were estimated by viscometric measurements. Enzyme activity was measured using cytidine 2', 3'-cyclic monophosphate as a substrate. The results from the measurements of distances between N(epsilon2) of His(12) and N(delta1) of His(119) and Stokes radius suggest (i) that the presence of sulfate ions decreases the distance between the catalytic His residues and increases the globular compactness, and (ii) that there is an expansion of the enzyme surface as well as elongation of the catalytic center in the presence of thiocyanate ion. These findings are in agreement with activity measurements

Temperature Dependence of Backbone Dynamics in Human Ileal Bile Acid-Binding Protein: Implications for the Mechanism of Ligand Binding

Human ileal bile acid-binding protein (I-BABP), a member of the family of intracellular lipid binding proteins plays a key role in the cellular trafficking and metabolic regulation of bile salts. The protein has two internal and, according to a recent study, an additional superficial binding site and binds di- and trihydroxy bile salts with positive cooperativity and a high degree of site-selectivity. Previously, in the apo form, we have identified an extensive network of conformational fluctuations on the millisecond time scale, which cease upon ligation. Additionally, ligand binding at room temperature was found to be accompanied by a slight rigidification of picosecond-nanosecond (ps-ns) backbone flexibility. In the current study, temperature-dependent N-15 NMR spin relaxation measurements were used to gain more insight into the role of dynamics in human I-BABP-bile salt recognition. According to our analysis, residues sensing a conformational exchange in the apo state can be grouped into two clusters with slightly different exchange rates. The entropy-enthalpy compensation observed for both clusters suggests a disorder-order transition between a ground and a sparsely populated higher energy state in the absence of ligands. Analysis of the faster, ps-ns motion of N-15-H-1 bond vectors indicates an unusual nonlinear temperature-dependence for both ligation states. Intriguingly, while bile salt binding results in a more uniform response to temperature change throughout the protein, the temperature derivative of the generalized order parameter shows different responses to temperature increase for the two forms of the protein in the investigated temperature range. Analysis of both slow and fast motions in human I-BABP indicates largely different energy landscapes for the apo and halo states suggesting that optimization of binding interactions might be achieved by altering the dynamic behavior of specific segments in the protein

Neuer Kopf, alte Ideen? : "Normalisierung" des Front National unter Marine Le Pen

In this article, it is investigated whether vibrational entropy (VE) is an important contribution to the free energy of globular proteins at ambient conditions. VE represents the major configurational-entropy contribution of these proteins. By definition, it is an average of the configurational entropies of the protein within single minima of the energy landscape, weighted by their occupation probabilities. Its large part originates from thermal motion of flexible torsion angles giving rise to the finite peak widths observed in torsion angle distributions. While VE may affect the equilibrium properties of proteins, it is usually neglected in numerical calculations as its consideration is difficult. Moreover, it is sometimes believed that all well-packed conformations of a globular protein have similar VE anyway. Here, we measure explicitly the VE for six different conformations from simulation data of a test protein. Estimates are obtained using the quasi-harmonic approximation for three coordinate sets, Cartesian, bond-angle-torsion (BAT), and a new set termed rotamer-degeneracy lifted BAT coordinates by us. The new set gives improved estimates as it overcomes a known shortcoming of the quasi-harmonic approximation caused by multiply populated rotamer states, and it may serve for VE estimation of macromolecules in a very general context. The obtained VE values depend considerably on the type of coordinates used. However, for all coordinate sets we find large entropy differences between the conformations, of the order of the overall stability of the protein. This result may have important implications on the choice of free energy expressions used in software for protein structure prediction, protein design, and NMR refinement

A Rigidifying Salt-Bridge Favors the Activity of Thermophilic Enzyme at High Temperatures at the Expense of Low-Temperature Activity

Although enzymes from thermophiles thriving in hot habitats are more stable than their mesophilic homologs, they are often less active at low temperatures. One theory suggests that extra stabilizing interactions found in thermophilic enzymes may increase their rigidity and decrease enzymatic activity at lower temperatures. We used acylphosphatase as a model to study how flexibility affects enzymatic activity. This enzyme has a unique structural feature in that an invariant arginine residue, which takes part in catalysis, is restrained by a salt-bridge in the thermophilic homologs but not in its mesophilic homologs. Here, we demonstrate the trade-offs between flexibility and enzymatic activity by disrupting the salt-bridge in a thermophilic acylphosphatase and introducing it in the mesophilic human homolog. Our results suggest that the salt-bridge is a structural adaptation for thermophilic acylphosphatases as it entropically favors enzymatic activity at high temperatures by restricting the flexibility of the active-site residue. However, at low temperatures the salt-bridge reduces the enzymatic activity because of a steeper temperature-dependency of activity

Structure and Behavior of Human α-Thrombin upon Ligand Recognition: Thermodynamic and Molecular Dynamics Studies

Thrombin is a serine proteinase that plays a fundamental role in coagulation. In this study, we address the effects of ligand site recognition by alpha-thrombin on conformation and energetics in solution. Active site occupation induces large changes in secondary structure content in thrombin as shown by circular dichroism. Thrombin-D-Phe-Pro-Arg-chloromethyl ketone (PPACK) exhibits enhanced equilibrium and kinetic stability compared to free thrombin, whose difference is rooted in the unfolding step. Small-angle X-ray scattering (SAXS) measurements in solution reveal an overall similarity in the molecular envelope of thrombin and thrombin-PPACK, which differs from the crystal structure of thrombin. Molecular dynamics simulations performed with thrombin lead to different conformations than the one observed in the crystal structure. These data shed light on the diversity of thrombin conformers not previously observed in crystal structures with distinguished catalytic and conformational behaviors, which might have direct implications on novel strategies to design direct thrombin inhibitors

Peptide Conformer Acidity Analysis of Protein Flexibility Monitored by Hydrogen Exchange†

ABSTRACT: The amide hydrogens that are exposed to solvent in the high-resolution X-ray structures of ubiquitin, FK506-binding protein, chymotrypsin inhibitor 2, and rubredoxin span a billion-fold range in hydroxide-catalyzed exchange rates which are predictable by continuum dielectric methods. To facilitate analysis of transiently accessible amides, the hydroxide-catalyzed rate constants for every backbone amide of ubiquitin were determined under near physiological conditions. With the previously reported NMR-restrained molecular dynamics ensembles of ubiquitin (PDB codes 2NR2 and 2K39) used as representations of the Boltzmann-weighted conformational distribution, nearly all of the exchange rates for the highly exposed amides were more accurately predicted than by use of the high-resolution X-ray structure. More strikingly, predictions for the amide hydrogens of the NMR relaxation-restrained ensemble that become exposed to solvent in more than one but less than half of the 144 protein conformations in this ensemble were almost as accurate. In marked contrast, the exchange rates for many of the analogous amides in the residual dipolar coupling-restrained ubiquitin ensemble are substantially overestimated, as was particularly evident for the Ile 44 to Lys 48 segment which constitutes the primary interaction site for the proteasome targeting enzymes involved in polyubiquitylation. For both ensembles, “excited state ” conformers in this active site region having markedly elevated peptide acidities are represented at a population level that is 102 to 103 abov